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Magnetic inelastic neutron scattering studies of iron-based superconductors reveal a strongly 
temperature-dependent spin-fluctuation spectrum in the normal conducting state, which develops 
a prominent low-energy resonance feature when entering the superconducting state. Angle-resolved 
photoemission spectroscopy (ARPES) and scanning tunneling spectroscopy (STS) allow us to study 
the fingerprints of fluctuation modes via their interactions with electronic quasiparticles. We cal- 
culate such fingerprints in 122 iron pnictides using an experimentally motivated spin- fluctuation 
spectrum and make a number of predictions that can be tested in ARPES and STS experiments. 
This includes discussions of the quasiparticle scattering rate and the superconducting order parame- 
ter. In quantitative agreement with experiment we reproduce the quasiparticle dispersions obtained 
from momentum distribution curves as well as energy distribution curves. We discuss the rele- 
vance of the coupling between spin fluctuations and electronic excitations for the superconducting 
mechanism. 
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I. INTRODUCTION 

Shortly after the discovery of high-temperature super- 
conductivity in Fe-based pnictide compounds by Hosono 
and collaborator^ a magnetic Cooper-pairing mecha- 
nism was proposed.!^ A pure electron-phonon interac- 
tion as the primar y pa iring mechanism was, however, 
found to be unlikely!^ It is believed that the proximity 
of antiferromagnetic spin density waves and supercon- 
ductivity in the phase diagrarrP^ is likely to support a 
pairing scenario where superconductivity is driven dom- 
inantly by an electron-electron interaction mediated by 
spin fluctuations.^ ^ The spin excitation spectrum of the 
iron pnictides (and iron chalcogenides) shows pronounced 
similarities with other superconductors where a magnetic 
pairing mechanism is under debate. This includes for in- 
stance a strong temperature dependence of the spectrum 
in the normal state, which was studied in iron pnictides 
by Inosov et al}^ ^ as well as the presence of a spin res- 
onance feature below a gapped continuum in the super- 
conducting st ate, I^^HI^ which also appears in cuprates and 
in some heavy fermion superconductors. 

Angle-resolved photoemission spectroscopy (ARPES) 
experiments on iron pnictides^^ reveal a sharp Fermi sur- 
face consisting of multiple electron and hole pockets (see 
e.g. Ref. .14. and references herein). These exhibit com- 
parable superconducting order parameter amplitudes at 
Fermi surface sheets which are nearl y nested by the an- 
tiferromagnetic wave vector, Q .115116] It thus seems likely 
that magnetic and electronic order are closely related.^ 
In the case of antiferromagnetic spin fluctuations as the 
origin for pairing, the order parameter must have differ- 
ent sign at the electron and hole pockets. In this 
context, we recently endeavored to answer the question 



of how strongly electrons couple to spin fluctuations in 
these compounds.^ To this end, we investigated low en- 
ergy dispersion anomalies, whose position and shape can 
be traced back to a coupling of bosonic modes - a method 
which has been proven to be a powerful tool previously 
in the case of cuprate superconductors.^^ In iron-based 
superconductors such anomalies have been observed in 
the hole-doped 122 compound Bai_a^Ka^Fe2As2 by Wray 
et alW, Koitzsch et a/.,l^and Richard et a/.l^ These ex- 
periments reveal that the electronic dispersion features a 
self-energy effect, which is most pronounced at an energy 
eo ~ 25 meV. 

Shortly before that, inelastic neutron scattering stud- 
ies by Christianson et alM' of the same compound re- 
vealed the development of a spin resonance in the su- 
perconducting state. This resonance appears at an en- 
ergy Qres ^ 14meV, is situated at the in-plane antifer- 
romagnetic wavevector Q = (tt, tt) and is weakly disper- 
sive in the Qz direction. ^^ ^^ Guided by similar studies in 
cuprate superconductors, we look for corresponding sig- 
natures in the available ARPES data.l^^H^ There, it was 
observed that (i) the superconducting excitation gap at 
the Fermi surface pockets nested by Q has the absolute 
value A ^ 12meV, (ii) the self-energy effect occurs at 
an energy eo ^ ^res + A, and (iii) eo follows an order- 
parameter-like evolution in temperature as the resonance 
energy does. Thus, it is intuitively plausible to conjec- 
ture that the magnetic resonance is responsible for the 
observed self-energy effects. By numerical calculation we 
show that this conjecture is theoretically well founded. 

Motivated by the experimental procedures'^ to quan- 
tify self-energy effects in angle-resolved photoemission 
experiments, we concentrated in a previous Letter on the 
so-called effective self- energy as extracted directly from 
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the electronic spectral function!^ Such an approach al- 
lowed us to make an immediate comparison with ARPES 
experiments. In the present work we are taking a more 
detailed look into our model, introducing the formalism 
used for our calculations, discussing the renormalization 
effects entering the electronic structure as well as the 
Fermi surface, and applying our model to photoemission 
and tunneling experiments. 



II. THE SPIN-FLUCTUATION SPECTRUM 
AND THE RESONANCE MODE 

The appearance of a low-energy resonance mode in the 
dynamic magnetic susceptibility upon entering the super- 
conducting state is a feature well known from unconven- 
tional superconductors (for a review concerning cuprates 
see Ref. 21). As it results from scattering between 
Fermi surface sheets with opposite order-parameter sign 
its detailed experimental and theoretical investigation 
is a powerful tool to s tudy the symmetry of the su- 
perconducting gapj -^^l^^ l Various microscopic theoretical 
treatments are able to reproduce the resonance feature, 
such as, for example, in methods employing the ran- 
dom phase approximation (RPA) or in the fluctuation 
exchange (FLEX) approximation. However in all known 
treatments specific assumptions and approximations are 
necessary. Strong correlations usually lead to vertex cor- 
rections in a Feynman diagrammatic technique, which, in 
most cases, are either neglected or, in the absence of any 
good expansion parameter, introduced by ad hoc meth- 
ods or educated guesses (as in conserving approximations 
such as FLEX). As our results do not depend on the 
specific microscopic model for the resonance excitation, 
we prefer here a semiphenomenological approach. Most 
of the magnetic correlations are included automatically 
within the semiphenomenological approach used here by 
relying on experimental results. 

In the following sections we introduce an effective 
model for the spin susceptibility and establish the con- 
nection to recent neutron scattering experiments. 



A. Spin susceptibility 

In this section we summarize the semiphenomenologi- 
cal approach we employ in this paper, which has already 
been successfully used to describe optim ally doped and 
overdoped cuprate superconductors.^^ 29 32 motivate 
its application to iron-based superconductors we recall 
one possible interpretation of the resonance mode in the 
spin susceptibility. We underline, however, that this in- 
terpretation is not necessary for our predictions to hold, 
as we base the spin excitation spectrum on the experi- 
mental observations. 

Let us divide the phase space of electronic excitations 
into high-energy and low-energy regions, the latter being 
located around the Fermi surfaces and populated by the 



electronic quasiparticles (see e.g. Ref. |33j). Incoherent 
spin fluctuations are dominated by high-energy electronic 
excitations; however, fine features such as the resonance 
mode result from modifications in the low-energy elec- 
tronic spectrum. In the case of Fermi surface nesting in 
the vicinity of certain "nesting points" a common approx- 
imation is given by the RPA enhanced susceptibility 



{Xhigh(9) 



h(9)}"\(i) 



where rj^.g^^(gr) are vertex functions involving only 

high-energy processes, II(a;, q) describes the polarization 
due to low-energy processes, and a hat denotes a ma- 
trix structure due to orbital degrees of freedom. (note 
that the vertex functions occur twice here as they con- 
tain no low-energy inclusions). We neglect for simplicity 
in the following a possible gr-dependence of the vertex 
functions. The high-energy part Xhigh accounts for inter- 
mediate or long-range antiferromagnetic correlations and 
thus can be well approximated by the phenomenological 
Ornstein-Zernike form Xhigh(g') = Xq / {I ^ ^'^ \q - Q\'^) 
with ^ = ^(Q) a correlation length, and where xq is 
nonzero only for the relevant orbitals. 

In order to find a suitable analytic form that fits the ex- 
perimental data near the antiferromagnetic wave vector 
Q, it is a common procedure to start from an approxima- 
tion for Il{uj^Q) in the energy region of interest. To this 
end, let us consider a simple two band model consisting 
of holelike and electronlike bands, ("h and respectively, 
nearly nested by the antiferromagnetic wave vector Q, 
and let us neglect for a moment the orbital structure. 
In the normal state interband scattering can lead to a 
relaxation process exciting a particle-hole pair around 
the Fermi surface. Such excitations most likely ap- 
pear near the nesting points at the electron (e) and 
hole (h) pockets, i.e. where Ch(^) ^ Ce{k — Q) -\- uo. 
Linearizing the dispersions around the nesting points, 
[Ch(^),Ce(^-Q)] = [kj:,ky\V^ lead to lmIi{uj,Q) ex 

rShrSh / ^Ch / <e det( v-1) [/(Ch) - /(Ce)] - [Ch - 

Ce]), where / is the Fermi distribution function, and V 
is a matrix with components proportional to the Fermi 
velocity components. For small frequencies ^ the 
imaginary part is linear in energy (and the real part neg- 
ligibly small), leading to an approximation for the sus- 
ceptibility 



X'(^,9) 



XQ 



(2) 



A detailed investigation of the spin dynamics in op- 
timally doped BaFei.85Coo.i5As2 was performed by In- 
osov et al^^ By including a temperature dependence in 
the parameters in Eq. |2|, that is replacing xq ~^ Xt-, 
and l^max ^ ^^x? ^^^Y were able to show that 
the energy and momentum dependence as well as the 
temperature behavior of Eq. ([2| fit well with the normal 
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state behavior, which vahdates theoretical models based 
on an itinerant description of magnetic excitations. 

In the superconducting state a low energy resonance 
appears with a weight that follows the temperature de- 
pendence of the superconducting gap. This again can be 
understood in a similar manner. Particle-hole excitations 
in the superconducting state that result from scattering 
between the two bands are suppressed below the two- 
particle excitation gap 



optimal doping the resonance energy is centered around 
the antiferromagnetic wave vector, whereas at higher 
dopings the spectral weight becomes large around two 
wave vectors and is therefore incommensurate.^^ Within 
the accuracy of the known parameters, neglecting the 
dispersive features will be sufficient for our purpose. 



B. Bai-xKxFesAss 



2Aq = min|fceF5} {\^h{k)\ + \A,{k - q)\) , (3) 

where Ae/h is the superconducting gap at the 
electron/hole pocket. Excitations in the contin- 
uum set in above this threshold, i.e., Imn(cj,gr) ex 
2Aq sign(a;) — 2Aq). For different signs between 
Ah (A;) and Ae(/c — g), the coherence factors appearing 
in 11(0;, q) are such that the real part has a logarithmic 
divergence at 2Aq. Furthermore, it can be expanded for 
small cj, i.e., Ren(cj, g) ex uj'^/2Aq. We define the reso- 
nance energy n.^s^q via Xhigh(<?)rgigh n(a;, q) ^ a;V^L,q 
and insert this expression into Eq. ([T]) to get 



Xq 



(4) 



res,q 



Here we have introduced a small broadening Fres <^ 
^res,q of the resouaucc mode, which accounts for its ex- 
perimentally observed finite width. The resonance mode 
can be identified with the pole in Eq. Q, which is 
moved by the amount Fresinto the lower complex half 
plane. If the correlation length is sufficiently large, 
the susceptibility is sharply peaked around Q, in the re- 
gion 

We rewrite Eq. Q as 



Wq/tT 



2Qyqs{q) 



r,(5) 



with Qresiq) = ^res,q(l + ^tIQ " Q?)^^^ and Wq = 

7rxQ^res,Q/2. In the following we adopt the approxi- 
mation to neglect the dispersion of the resonance energy 
with momentum. This is, we replace Eq. ([5| by 



Wq/tT 



{UJ + zFres) 



2 • 



(6) 



This form has the advantage of fitting the energy and 
momentum width independently. From experiment is it 
known that the dispersive features have much less weight 
then the resonance itself, and thus this will lead to a 
strongly reduced infiuence on the electronic single parti- 
cle dispersions. There have been experimental and the- 
oretical investigati ons o f the doping dependence of the 
bosonic dispersion.'^^'^ They show that the incommen- 
surability of the resonance and the dispersion of the reso- 
nance mode depend on the doping level. Below and up to 



In this section we outline our phenomenological ap- 
proach for the case of iron based superconductors. We 
shall focus on hole-doped Bai_xKxFe2As2, since we will 
compare our predictions with the available experimental 
angle-resolved photoemission spectra for this compound. 

Spin fiuctuations that originate from the spin of the 
conduction electrons are closely connected to the Fermi 
surface topology of the electronic structure. In our model 
we employ a tight-binding fit in an orbital basis that was 
obtained from the density functional theory (DFT) band 
structure of BaFe2As2 by Graser et al*^ The dominant 
contribution to the density of states in the energy range 
±1-2 eV around the chemical potential originates from 
the five Fe 3d orbitals (with some hybridization with the 
As 4p orbitals mainly for energies above the chemical 
potential). This energy range will be sufficient for our 
purposes. The resulting reduction of the Hamiltonian to 
a five-orbital basis reads 

Smn^n]dn{ka). (7) 

ka mn 

Here d\^(k(j) creates an electron with momentum k and 
spin (J in the orbital m, where m = 1, 2 • • • 5 corresponds 
to the five orbitals dxz^ dyz^ dj;2_y2^ dxy^ d^z2_j,2. The 
parameters ^mn and are listed in Ref. 36, The canon- 
ical transformation 



diika) = ^ a^{k)dl{ka) 



(8) 



diagonalizes the Hamiltonian, leading to eigenvalues 
C^{k) and eigenvectors a^ik) = {mk\iik)^ where ji rep- 
resents the band index. Because Cmn{—k) = Cmn(^)* the 
eigenvectors can be choosen such that a^{—k) = a^ik)* 
holds. The set of eigenvectors for each k is ort honor mal 
(or can be chosen so), i.e. Xlm ~ and 

Ep«I^(fc)*«;^(fc) = '5„^n. 

In order to simulate hole-doping we apply a rigid shift 
of the chemical potential, i.e. we define the new band 
structure C/i(^) = C^(^) ~ Mo with jiQ = — 50meV. This 
was shown to be at least applicable for materials in which 
doping occurs via substitution in regions outside of the 
conducting Fe-As planes and subsequent charge trans- 
fer, as in the case of Bai_xKxFe2As2.^^ This assump- 
tion is consistent with our results for weak to moder- 
ate coupling, which show that the renormalization of the 
chemical potential due to the coupling between electronic 
excitations and spin fiuctuations is in this case well ap- 
proximated by a linear dependence between the chemical 
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potential and the charge carrier concentration (for details 
see below). 

The dominant processes for spin fluctuations are the 
interactions between electrons at the various Fe orbitals 
of the same atom (inter- and intra-orbital Coulomb inter- 
action, Hund's coupling, and intra-orbital pair hopping). 
Thus, spin fluctuations are described naturally in an or- 
bital basis, whereas electronic excitations are easier to 
discuss within a band picture. Because we are interested 
in the influence of spin fluctuations on electronic excita- 
tions, it is useful to recall the orbital characters at each 
Fermi surface point. 

Throughout this paper, we adopt the notation of 
Ref. [36l corresponding to 1 Fe/unit cell. Fig. [l] summa- 
rizes the orbital contributions to the Fermi surface for 
three values of /c^: 0, 7r/2, and tt. 

Starting with kz = 0, Figs. [TJa)-[TJc) show the partial 
spectral function for the bare bands. 



(a) 



An{k) 



E 



(9) 



for the orbitals n = dy^^ dj.2_y2 and dxy^ respectively 
(where 5 = 4meV is an broadening factor), whereas 
Fig. [TJd) shows the dominant orbital contributions to 
the Fermi surface only. Corresponding orbital characters 
for A;;^ = 7r/2 and kz = tt are shown in Figs. [TJe)-jT][h) and 
[lji)-[ljl), respectively. 

Magnetic excitations can be described in terms of the 
dynamic spin susceptibility x^^? where the measured 
susceptibility x is the sum over all orbital contribu- 
tions, X — X^mn/^mT- ^au be seen from Fig. [l] the 
hole pockets at (0,0,/^;^) and (7r,7r, /c^) are nearly nested 
(taking into account the finite momentum width of 
the spin excitation) to the electron pockets at (0,7r, A:^) 
and (7r,0,/C;^) by the wave vectors Q\ = (ii^^^qz) and 
Q2 = (O^Tr,^^), which connect Fermi surface sheets with 
same orbital character and correspond to antiferromag- 
netic correlations in the respective direction. Fluctuation 
exchange approaches have shown that the magnetic mode 
predominantly scatters between orbitals with the same 
character. ^^ ^^ According to that we neglect all inter- 
orbital contributions in our model, i.e. we only couple 
the part Xn = Xnn electronic excitations, and we 

assume X « EnXn- 

As can be seen from Fig. [ij the dx2_y2 orbital has 
negligible intra-orbital contributions to nesting, because 
the electron pockets have very little of its character. 
Thus, we only take into account the remaining orbitals 
(n = dxz^ dyz, dxy = 1, 2, 4). To simplify notation, we 
introduce the parameters 



1/2 for(na) = (21),(12),(41),(42) 
else 



(10) 



to account for the orbital selective mode coupling via the 
wave vectors <5a=i,2 (see Fig. [ijl)). With this, the nor- 
mal state susceptibihty will be modelled by an equivalent 



(c) 




(b) 



(e) 



(g) " 



dyz 


1 




^ 




r 


> 




1 













K 




dxy 










5 












n 





7. 




K 





250 
200 
150 




250 
200 



(k) 


r 

dxy 






^0 


D 





C 















K 


7 



250 
200 



dx-^ —y^ 



(d) 




(f) 



(h) 




(j) " 



(I) 




FIG. 1: Orbital character at the Fermi surface for 
Bai-ccKxFe2As2 (x = 0.4): (a)-(c) Partial spectral function 
An{k) = E;. K{k)\^S/{6^ + C^{k)^) at k, = for the orbitals 
n = dyz, dx2_y2 and dxy, respectively. The corresponding fig- 
ure for the dxz orbital is obtained from (a) by a 90"^ rotation, 
(d) Dominant orbital contributions to the Fermi surface at 
kz = (with the dxz orbital in green). The various Fermi 
surface sheets are named a 1,2 (hole sheets), /3i,2 (electron 
sheets), and 7 (hole sheet), (e)-(h) the same as in (a)-(d) for 
kz = 7t/2. (i)-(l) The same as in (a)-(d) for kz — ti. The two 
nesting vectors Qi and Q2 are shown in (1). 
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of the Ornstein-Zernike form in Eq. ([2|, 

a=l,2 



, (11) 



with temperature-dependent parameters Xt — Xo/{T -\- 
= ^o/VT^O, and Q^^^^ = Qo{T^O), as motivated 
by Ref. ^1 Here is the Curie- Weiss temperature. 

In the superconducting state a particle-hole excitation 
gap opens up in the spin excitation spectrum, a reso- 
nance peak appears within this gap, and spectral weight 
is shifted into this magnetic resonance. Accordingly, the 
magnetic spectrum now consists of two parts, the low- 
energy resonance and the particle-hole continuum. 

Particle-hole excitations appear above a temperature 
dependent threshold of 



2A(T) 



|A^(/c^)|), 



(12) 

where and are the superconducting gaps at the 
nested pockets, and and are the corresponding 
Fermi wave vectors. For high excitation energies the sus- 
ceptibility should recover the normal state behavior. Mo- 
tivated by these observations we approximate the gapped 
continuum by the same functional form (but with dif- 
ferent magnitude; see below) as in the normal state, 
Eq. ( 11 ), and write for the imaginary part of the complex 



dynamical susceptibility 



for 



> 2A(T). (13) 



Neutron scattering experiments have shown that the res- 
onance follows an order parameter like evolution. There- 
fore, we assume the temperature dependence 



(14) 



and we apply the Lorentzian form of Eq. ([5| in order to 
model the resonance below the threshold, i.e. 



xT{^.Q)=x7{^.Q) for |c^|<2A(T), (15) 



with 

Xni^^Q) = Yl 



a=l,2 



TT 



^^res 



(16) 



Here the momentum distribution enters via the weight 
function 



l + ^?\q-Qa\ 



(17) 



Note that the momentum and orbital dependence of the 
susceptibility is contained in the factors bn^a (Eq. 
We determine the real part (up to a constant, see be 
from x"^n t>y exploiting Kramers-Kronig relations 



10). 



ow 



1 



duj' 



//SC/ / \ 



UJ' 



■ UJ 



(18) 



The resulting energy dependence of the dynamical sus- 
ceptibility and its variation with temperature are shown 
in Fig. (2] 

The resonance weight factors in Eq. (pTj ) are 
determined by a local sum rule (Appendix |B|) which 
fixes the ratio between resonance and continuum, 
w'^/Xt- They are chosen so that the integrated spin 
structure factor, J^^du J dqS{uo^q)^ with S{uo^q) = 

2^Xln Xn{^-> 9)/(l ~ e"^'^/^^-^), remains temperature in- 
dependent. In general the weight factors Xt 
pend on the orbital character. In order to determine the 



relative weight between 



{x¥ = Xt ) 



(Xt^) we studied the electronic dispersions at the X point 
of the Brillouin zone and find that experiments are best 
reproduced when the weight is roughly equal. We thus 



set w^^ = = w> 
and are left with 



xy 



wt and Xr 



yz 

Xt 



xy 

Xt 



Xt 



2A(T) 



duo 



2A(T) 



d\ 



X"n{^^Q) - X"1{^^Q) 
1 _ p-hoj/kBT 



= (19) 



to fix wt/xt- We illustrate this procedure in Figs. [2|a) 
and [2|c): When integrating over energy and momen- 
tum the green area above and under the dashed curve 
has to be (approximately) the same. The ratio is in 
excellent agreement with the functional form wt/xt = 
{l-TyT^)wT=0K/XT=0K (Fig-lJe)). 

We apply the model above to Bao.6Ko.4Fe2As2 and use 
the experimentally motivated parameter set in Table [T[ 
The only remaining parameter to be fixed is the quan- 
tity g'^Xo^ where the overall weight xo can in principle 
be extracted from experiment but is not known for this 
compound, and g is the coupling constant between elec- 
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FIG. 2: Energy dependence of the dynamic spin susceptibil- 
ity in units of xo (lO'^K)"^: (a) and (c) illustration of the 
local sum-rule at two different temperatures, (b) Tempera- 
ture evolution of the magnetic resonance, (c) Temperature 
dependence of the resonance weight factor wt / Xt • 
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TABLE I: Parameter set appropriate for Bao.6Ko.4Fe2As2 



^^0 


Co 


e 


& 




-L res 


0.375 


5.84 K^/^ 


30 K 


2 


15.5 meV 


3 meV 



trons and spin fluctuations. To account for the periodic- 
ity in momentum space we replace the factors \q — Q^P 
in Eq. ^ and (fTTl) by 



Qc 



4 



sm 



sm 



Qy Qc 



(20) 

and ne glect the g'^-dependence for our purpose as it varies 
weaklyPMl We use the in-plane lattice constant a and 
the out-of-plane lattice constant c as units of length. 

The resulting energy dependence of the spin-excitation 
spectrum and its evolution with temperature is presented 
in Fig. [2] (b). Experimentally the resonance appears at 
an energy of ^^7^^ ^ 14meV'I2 below the particle-hole 
continuum u < 2A(T = 15 K) ^ 24meV.^^ In momen- 
tum space the mode is peaked around the wave vectors 
Qa=i,2 with a correlation length of nearly twice the 
lattice constant. When the temperature is decreased be- 
low Tc, the resonance gains weight and eventually be- 
comes the dominating part of the spectrum at low ener- 
gies. Note that the gain in weight is due to two parallel 
processes: a suppression of high-energy weight due to an 
increase in correlations (already present in the normal 
state), and an increase of the superconducting gap (and 
thus the spin excitation gap) with decreasing tempera- 
ture. 

As we will see later, the characteristic temperature de- 
pendence of the resonance mode imprints itself onto the 
electronic spectrum, where the resonance leads to a sig- 
nificant effect below 30 K. 



A. Formalism 

In this section we summarize the formalism we use 
to study these effects. We use a perturbative approach 
based on Gor'kov Green's functions with the coupling 
between the conduction electrons and the spin fiuctu- 
ations as the expansion parameter. The unperturbed 
Green's functions are diagonal in band index, with nor- 
mal (diagonal) and anomalous (off-diagonal) components 

G^^^ (e, k) and F^^^ (e, k). The renormalized Green's func- 
tions, Gjj^u and F^y are not diagonal in band index due 
to the interactions induced by spin-fiuct nation exchange. 
The Green's functions in the orbital basis, Gmn and Fmn^ 
are related to those in band representation, Gjj^u and F^^y, 

by 

Gmn{e,k) = Y,<ik)<{krG^Ae,k), (21) 
Fmn{e,k) = ^a';{k)a:;{-k)F^,{e,k). (22) 

We couple electrons to the spin-fiuctuation spectrum 
with an energy and momentum independent coupling 
constant g (instantaneous and local coupling). The re- 
tarded diagonal and off-diagonal self-energies are then 
given by 



R 

mn 



(23) 



written in terms of the retarded (R) and Keldysh (K) 
Green's functions 



9^ (^nn 



^ nn An / ' 



tl V nn An ^ nn An / ' 



(24) 
(25) 



with {A * B){e, fe) = / / (|^A(e -LO,k-q) B{co, q) 
as explained in Ref. 21, In equilibrium the Keldysh com- 
ponents are given by 



III. COUPLING TO SPIN FLUCTUATIONS 

We are ultimately interested in the renormalization 
of the electronic dispersion and in the electronic life- 
time as a result of the coupling of conduction electrons 
to spin fiuctuations both in the normal and in the su- 
perconducting state. We model this coupling by an ef- 
fective electron-electron interaction mediated by the ex- 
change of spin fiuctuations. We concentrate on the lead- 
ing (quadratic) contribution in the coupling constant. 
This approximation can be justified in part by a phase- 
space consideration. Small parameters are introduced by 
the restricted phase-space areas available for electronic 
quasiparticle excitations as well as for the low-lying spin- 
fiuctuation modes. This introduces stringent kinematic 
restrictions.!^ 



Xn(^,9) 

G^^(e, k) 
FL{e,k) 



2ilmx^{co,q) [l + 2 5(w)], 
2*ImGl(e,fc) [l-2/(e)], 
[F„f^„(e,fc) - F„f^„(-e, -fe)] [1 - 2/(e)], 



where / and b are the fermionic and bosonic distribu- 
tion functions. The normal and superconducting state 
dispersion relations are given by Xn = Xn Xn = Xn 
respectively, where n = 1,2, 4. The self-energies then 
enter the Dyson equation in terms of Nambu-Gor'kov 
Green's functions 



A(0)R 
mn 

E 



mn 



(26) 



where we used the completeness relation 
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and 



- f 



R 

mn 

yR 

^mn 

mn 



T^R 

mn 
_fiR 
^mn 

mn 



(27) 
(28) 



with A(e, k) = A*{-e, -k). From Eq. (|26| we obtain 

X [Z„(e, fe)(e + i5)i - C;(e, fe)f3 - A"(e, fe)fi] , 

where 1 is a 2 x 2 unit matrix and {i = 1,2,3) are 
the Pauh matrices in Nambu space. The renormahzed 
dispersion and order parameter as weh as the renormal- 
ization function are given by 



c;(e,fe) 

A"(e,fe) 
Z„{€,k) 



Uk)^?^i^^^^, (30) 



S^(e,fc)-S^(e,fc) 
2(e + id) 



(31) 
(32) 



using the fact that in our case = Note that the 
real part of the dynamical susceptibility is determined by 
a Kramers-Kronig analysis only up to a constant. This 
would lead to energy independent contributions to (J^ 
and A"^. However, any such contribution can be thought 
about as absorbed into the band structure and A^, which 
enter in our approach as phenomenological parameters. 
Finally, by inverting Eq. (29) numerically the spectral 



function is obtained from the first diagonal element of 
Eq. (I27l), 



A(6,/c) = --Im^G^^(6,/c). 



(33) 



For analytic properties of the self-energies and a sum rule 
we refer to Appendix [Dj The density of electrons is given 
in terms of the spectral function by 



2 
V 



de 



(2^ 



f{e)A{e,k), 



(34) 



When switching on the interactions, it must be ensured 
that the chemical potential is adjusted such that the elec- 
tronic density stays constant. 



We calculate the convolutions in Eqs. (24) and (25) nu- 



merically by a fast Fourier transform using bare Green's 
functions (a procedure supported by the numerical stud- 
ies in Ref. 43) with a broadening parameter of 6 = 4 meV. 
For this we use a very fine 512 x 512 x 8 /c-mesh and 128 
points in energy space. In addition a high-energy cut- 
off of ujc = 200 meV was introduced in the spectrum of 
spin excitations. The exact value of this is, however, 
less important, as any change in the high-energy part 
contributes only to an additional renormalization factor 
which we address next. 



High-energy cutoff and high-energy 
renormalization 



As already mentioned, we use a high-energy cutoff 
ujc = 200 meV in the spin excitation spectrum. The pre- 
cise value of this cutoff is, however, not essential. To see 
this, let us assume that we change the high-energy cutoff 
from the value uJc to cJc+Acj. Accordingly, the self-energy 
will have two terms, which we call the low-energy and the 
high-energy part, i.e. = + 1]™. For sufficiently 
large cutoff cJc, the high-energy part contributes at low 
energies (well below Uc) mostly to the real part and is 



linear in energy, i.e. 



vHE 



HE 

e,k 



-hn e {hn e R) 



(see Appendix [C|. We exploit this fact to define an 
energy independent high-energy renormalization factor 
= 1 - S'HE/e = 1 + 6„. 
The new Green's functions and self-energies corre- 
sponding to the cutoff UJc + Acj are determined by the 
assignments Q ^ Q/^n^. ^ A^/Z^^, and 
1 + (Zn - 1)/Z™ in the expression for [G'^^]"\ and 

[Gmn]~^ Zn^[Gmn]~^' Thcsc relations simplify fur- 
ther for equal coupling constants for all orbit als, in which 
case Z^^ = Zhe is independent of the orbital index n, 
and G^^ ^ G^^/Zhe. 

We treat Zhe as a free parameter, which is of order 
1 and modestly temperature dependent. This tempera- 
ture dependence might be at first view a bit surprising, 
as from Fermi liquid theory one is used to temperature 
variations of the high-energy vertices that are negligible. 
However, one should remember that in our case the sus- 
ceptibility is strongly temperature dependent in the nor- 
mal state, which includes temperature-dependent shifts 
of spectral weight between low and high-energy. This is 
manifestly non- Fermi-liquid behavior. 

We fix the high-energy renormalization factor above 
a certain reference temperature Tref (which we chose as 
50K) in the normal state Zj^g^ = 1, and we determine 
^HE lower temperatures so that the superconducting 
and the normal state dispersion merge for high energies, 
as observed in experiment j ^^ * ^^^ Our numerical solutions 
show that Z^g varies slowly 

in the temperature range T = 10 - 50 (Fig. |3] shows 
an example for g'^Xo = 1-17 x 10^/i|eV K). We underline 
that this additional high-energy renormalization applies 
only for energies well below the spin-fiuctuation cutoff 

We note that the thus determined Zhe is lower than 
one below a temperature of 50 K. In general, the renor- 
malization function due to the entire spin-fluctuation 
spectrum must be larger than one for zero energy. How- 
ever, this refers to the sum of our model spectrum and the 
correction due to the high-energy cutoff. Our model sus- 
ceptibility overestimates the high-energy contributions, 
so that a negative correction at high energies is in place, 
leading to negative values for 6^; the sum Z^ + 6^ at zero 
energy is, however, always larger than one. 
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FIG. 3: The high-energy renormahzation factor Zhe shows a 
hnear dependence in temperature. 



IV. SELF ENERGY EFFECTS 

In this section we discuss effects resulting from the self- 
energies, Eqs. (24) and (25). We start with the diagonal 
self-energies, which determine the quasiparticle scatter- 
ing rate and the quasiparticle band renormahzation, and 
proceed then with the off-diagonal self-energies, which 
renormalize the superconducting order parameter. 

All results in this section are for ^^Xo = 1.17x 10^/i|eV 

K. 



and is given in the normal state for the respective orbital 

^ — dxzi dyzi dxy by 



r„(e,fc) 



Im{E^(e,fc)} 
Re{Z„(e,fc)} ■ 



(35) 



A. Scattering rate and band renormalization 



The influence of the spin resonance mode on electronic 
quasiparticles is most pronounced in the vicinity of the 
Fermi surface. Therefore the scattering rate should ex- 
hibit characteristics of the Fermi surface topology. In 
Figs.[4];a) andgb) we show the kx and ky dependence 
of the low-energy dyz^xy scattering rate in the /c^ = 
plane and we show the quasiparticle renormalization fac- 
tor Re{Zn} in Figs. [I]; c) andgd). 

In the case of the dxz.yz orbitals [Fig. [4]^a) and|4|b)], 
both feature a two-hump structure in the middle of the 
Brillouin zone and a even larger peak at Kyz = (±7r,0) 
and Kxz = (0,±7r). As we demonstrate in Fig. [5| 
the hump structure is a clear signature of the respec- 
tive dxz,yz-o^^^^^^ contribution on the Fermi surface. In 
Figs. |5|b) andjsjd), we have shifted the partial spectral 
functions shown in Figs. [TJa) and [TJc) by the antiferro- 
magnetic wave vector Qi in the case of the dyz orbital 
and by Qi and Q2 in the case of the dxy orbital. It can 
be seen that the two-hump-structure in Fig. [sja) corre- 
sponds to the dyz contributions of the /3i pockets which 
lie opposite and well separated. On the other hand, the 
peak feature at Kyz = (±7r, 0) occurs due to the hole- 
like Fermi surface sheets at the a pockets where the dyz 
contributions enclose the F = (0,0) point. The scat- 



The imaginary part of the diagonal self-energy, 
i;^(e,/c), determines the scattering rate of electronic 
quasiparticles. In the presence of a renormalization fac- 
tor Z^(e, /c), Eq. (32), the scattering rate is renormalized 






FIG. 4: Momentum dependence of the low-energy scattering 
rate Tnie 0, k) in the dyz {n — yz) and the dxy (n — xy) 
channel at T = 31K [(a) and (b)], as well as the momentum 
dependence of the renormalization factor Zn(e 0, k) for the 
respective orbitals n [(c) and (d)]. 




FIG. 5: Intra-orbital scattering channels: The significant fea- 
tures in the dyz and dxy scattering rate, shown in (a) and (c) 
correspond to the different orbital contributions to the Fermi 
surface. We show in (b) and (d) the contributions of the dyz 
and the dxy orbitals to the Fermi surface like in Fig. [l] but 
shifted by the scattering vector Qi in the case of dyz and 
Qi and Q2 in the case of dxy. The scattering rate reflects 
the underlying structure of the Fermi surface shifted by the 
antiferromagnetic wave vectors. 
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FIG. 6: Scattering rate at T = 15 K: Energy-momentum 
dependence along the cut {ky G [— tt : tt), /c^c = A;^ = 0} in 
the dyz and dxy channel [(a) and (c)]. Energy dependence of 
the scattering rate at fixed momentum as indicated at the left 
hand side by the dotted blue and full black line, respectively 
[(b) and (d)]. Scattering sets in at the energy Qres- Peaks 
appear at the energies Qres + Aq and Qres + which are due 
to the coupling of the resonance to excitations at the Fermi 
surface and the van Hove singularities at the hole pockets 
centered at (7r,7r,z), respectively. 



tering rate is broadened compared to the Fermi surface 
contributions simply because the susceptibility is broad- 
ened in momentum space. Looking at the dxy orbital in 
Figs.jsjc) andjsjd) we essentially observe the same char- 
acteristics. The broader peak at the F point results from 
the (3 pockets and the sharp ones at the K^y = (±7r, 0) 
and K^y = (0, ±7r) points originate from the 7 Fermi 
surface sheets at (±7r, ±7r). 

In the superconducting state the opening of a particle- 
hole excitation gap leads to a suppression of the particle 
decay below e < 2A(T). However, the appearance of 
the resonance allows for scattering processes below the 
continuum threshold. As a consequence, inelastic scat- 
tering sets in above the resonance energy, i.e., e > l^^g, 
as depicted in Fig. [gI where the scattering rates Tdy^ and 
F^^^ are shown for T = 15 K. The largest contributions to 
scattering arise from states near the Fermi surfaces. This 
means that the resonance imprint is most significant at 
an energy of ftres + Aq, where for each considered Fermi 
surface wave vector, /c, Aq denotes the gap averaged 
over a region of diameter around each of the points 
k±Q 1/2 (which are close to another Fermi surface nested 
to the original one) [Fig. [6]^b)]. An additional peak ap- 
pears in the da^^y "Scattering rate at an energy of (^res + ^7 
which is due to the van Hove singularity at the 7 pockets 
[Fig. [6]^d)]. In the following, we will refer to Qres + Aq 
as the kink-energy. In Appendix [A] we present similar 
investigations of the diagonal self-energies in Eqs. (30) 
and p2|). 



B. The superconducting order parameter 

The origin of the pairing instability may well be re- 
lated to the spin-fluctuation continuum, as demonstrated 
by recent FLEX calculations.^^ Our model is restricted 
to the low-energy region in the spin excitation spectrum. 
Consequently, we have to include the high-energy inco- 
herent part, which considerably contributes to pairing, 
separately. The energy range of interest (|cc;| < ujc = 
200 meV) gives only a partial contribution (about 40- 
50%) to the value of the superconducting order param- 
eter. Thus, we must add an additional contribution A^ 
resulting from the high-energy incoherent spin fluctua- 
tions. 

According to theoretical calculations t he or der param- 
eter dominant ly has a 5^ symmetry^-^^^IEl]^ which fol- 
lows approximately the form 



A^ = Ao cos(A:a,) cos{ky). 



(36) 



For our model we have chosen as a first initial guess 
A^ = A J, which is independent of the orbital character. 
In the absence of detailed experimental data, we resort to 
the temperature dependence Ao(T) = Ao(0) y^l — T/Tc, 
which is true in the BCS-limit in the vicinity of Tc, but 
deviates from the true value at much lower tempera- 
tures. ARPES measurements resolve a superconducting 
gap of A(15K) ^ 12meV at the inner hole-like a-pocket 
as well as at the nested electron pockets.— We choose 
Ao(0) = IS.lmeV so that the renormalized gap matches 
the experimental one at the particular points in the Bril- 
louin zone. 

The low-energy spin fluctuations then lead to an 
orbital-dependent part, <l>n according to Eq. ( [25|. We 
plot these and the gap in Fig. [?[ From FigT^c) it 
can be seen that in the dxz,yz channel the coupling via 
the Q12 = (±7r,0), (0, ±7r) wave vector favors an or- 
der parameter with opposite sign on the a and the Pi ^2 
pockets. The corresponding analysis for the dxy chan- 
nel in Fig. W^d) shows that the presence of the holelike 
pockets at (tt, tt, z) strongly support pairing with oppo- 
site signs between the holelike 7 pockets and the elec- 
tronlike Pi ^2 pockets and with equal signs between the 
ai^2 and 7 pockets. Taken together, all three orbitals 



promote an order parameter of the form of Eq. (36). 



by 



According to Eq. (29), the renormalized gap is given 



A^ + ^^(e,fc) 
Zn{e,k) 



(37) 



In Fig. [t] we compare the high-energy contribution A J, 
the off-diagonal self-energy <l>^, and the renormalized gap 
given by Eq. (37). We find that the spin- fluctuation 
spectrum in the energy range |cc;| < 200meV contributes 
nearly 40%-50% to the observed order parameter. To 
see this, note that at the F-point Z^y^ equals approxi- 
mately 2 [see Fig. [4jc)], thus the low-energy contribution 
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FIG. 7: Momentum dependence of the superconducting gap 
at e = OeV and T = 15 K as calculated with the ansatz 
in Eq. (36) for the high-energy contributions: (a) s^-gap. 
(c) and (dj show the off-diagonal self-energies in the dyz and 
dxy channel, (b) compares and the renormalized order 
parameter An = (A^ + $^)/Zn along the cut {kx G [— tt : 
tt), = /c^ = 0}. (e) and (f) show the renormalized order 
parameter An in the dyz and dxy channel. 



^dy^ must be of the same order as to arrive at the 
^dyz shown in Figs. [T^b) andjrj^e). Furthermore the su- 
perconducting gap originating from the low-energy part 
of the spin-fluctuation spectrum supports an state. 
We performed additional calculations for an extended 5- 
wave contribution [co^ikx) + cosiky)] to the order pa- 
rameter ansatz in Eq. (36). According to our findings, a 



large extended 5-wave contribution would imply a nodal 
structur e in th e partial density of states (for its definition 
see Sec. VB4). A RPE^ ^ and quasiparticle interference 
(QPI) experiment J^^EU however show a fully gapped or- 
der parameter. For our purpose it is sufficient to use the 
ansatz in Eq. (36). 

In AppendixjA we also discuss the energy and momen- 
tum dependence of the off-diagonal selfenergy, Eq. (31). 



V. THE ELECTRONIC SPECTRAL FUNCTION 

Self-energies are measurable via their effect on the line 
shape of the spectral function and the associated disper- 
sive features. The spectral function can be measured by 



ARPES techniques. In this section we will apply pro- 
cedures to our theoretically obtained spectral function 
that routinely have been applied in ARPES experiments. 
This allows us to compare our results directly with ex- 
perimentally observed quantities, such as effective self- 
energies extracted from ARPES measurements. In addi- 
tion we will present results that can be compared with 
other spectroscopic experiments such as c-axis tunnel- 
ing through superconductor-insulator-normal metal and 
superconductor-insulator-superconductor junctions. 



A. Normal state dispersion, Fermi surfaces, and 
chemical potential 

Dispersions obtained by ARPES experiments usually 
differ from density functional calculations by a renor- 
malization factor of roughly 1.5-2.^^ We show that 
this renormalization can be well accounted for by the in- 
coherent high-energy spin-fiuct nation spectrum coupling 
to the electronic excitations. In Fig. [8] we compare the 
spectral function obtained from the bare electronic struc- 
ture with the one obtained after coupling to the spin- 
fluctuation continuum at T = 50 K. As can be seen, there 
is a pronounced energy broadening of the spectral func- 
tion at high energies, whereas it remains sharp in the 
vicinity of the Fermi surface. The high-energy broad- 
ening increases with excitation energy. Furthermore, the 
figure shows that the energy bands are strongly renormal- 
ized, in particular near the X and M points. By choosing 
the coupling constant g to be the same for all orbitals we 
are able to reproduce the experimen tally observed shal- 
low electron pocket near the X point .'^^'^ 

The Luttinger theorem requires that the volume of the 
Fermi surfaces must stay constant as a consequence of 
particle conservation. In Fig. [9] (a) the chemical potential 
dfi necessary to maintain a constant particle number is 
shown. By shifting the bands, ({k^jno) = C{^) ~ Mo ^ 
("(fe, /i) = C{^) ~ the density of electrons p, given by 
Eq. ([34|, is fixed. In the following, we denote by 
the spectral function calculated for chemical potential /i 
and for self-energies ^[A^]. For the unperturbed band 
structure /io = — 50meV as chosen in Sec. 
to simulate hole doping. 



II B in order 



To obtain the correct chemical potential we need to 
find /i such that = po with = So(/i)], where 

Eo(/i) is calculated from the unperturbed spectral func- 
tions, A^^\ with chemical potential /i. In a first step we 
calculate self energies Ho(Mo) = So[^m?]- We then, in a 
next step, calculate A^^ = A[/io, So(mo)] and from this 
pi = In order to obtain a first approximation to 

the chemical potential, we linearize its functional form as 
function of the charge carrier density, defining 



Pi - Po 
dp/ dp ' 



(38) 



The partial derivative is obtained by making a second 
calculation for a shifted potential p = po -\- Ap with 



11 




r X M r 



FIG. 8: Spectral function along the cut in the first Brillouin zone as indicated in the inset. Here, T = 50K (which corresponds 
to the normal state reference temperature we have used for analyzing the spectral functions in the superconducting states as 
discussed in the text), and g^Xo = 1.17 x 10^/i|eVK. The black wide-dashed lines present the bare band structure. Zero energy 
corresponds to the unrenormalized chemical potential (/io, for the bare band structure), and the shift of the chemical potential, 
for the renormalized band structure is indicated as a second dashed line. 



A/i = 5 meV and repeating the procedure above leading to pi = p[Ajj] and 



(a) 0.06 
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FIG. 9: (a) Dependence of the chemical potential 6 /a on 
the coupling constant g^Xo- Shown is the result after the 
first, second, third, and fourth iteration at T = 50 K. (b) The 
electronic density of states for different coupling parameters 
= 0.08 X lOVIeVK and 1.17 x 10^1 eFK at T = 
(c) Fermi surfaces in a quarter of the Brillouin zone at 
for g^Xo = and 6/j. = and for g^Xo — 1-17 x 
lOVs eV K and 6fi = 22 meV, at T = 50 K. 



9^X0 
50 K. 



dp _ Pl-Pi 



(39) 



(0)i 

Ml J 



dp/dfi 



^^/ii = 

and so 



We then repeat these steps: Ho(mi) = ^o[A[ 

i;o(/il)] ^ p2= ^ P2= Pl 

forth until convergence. 

In Fig. [9Fa) we show Stij = pj — po for j = 1, 2, 3 as a 
function of g'^Xo- In Fig-H^b) we compare the density of 
states for two different coupling strengths for T = 50 K. 

Finally, in Fig. ^c) we show the Fermi surface for 



9 Xo = and for g^o — 1.17 x p^ eVK as obtained 
from the above procedure. The main effect is that all 
Fermi surface sheets shrink with increasing interaction, 
whereby a redistribution of populated states from the 
electronlike band to the holelike bands takes place near 
the Fermi surfaces. This has the most pronounced effect 
on the small hole pockets at the M points. We specu- 
late that for sufficiently strong spin fluctuations a Lif- 
shitz transition is triggered with an interaction-induced 
collapse of the entire hole pocket. 

Experimentally, it was shown that a shift of the 
chemic al po tential corresponds to a linear change in 
doping meaning that in the weak-coupling regime 
the coupling strength is linear to the doping level. This 
observation would correspond in our calculations to solu- 
tions that are already to a good approximation reached 
after the first iteration, pi ^ p. Our results sup- 
port this picture for the weak-coupling regime, g'^Xo < 
1 X lO^/i^eFK; however, they show pronounced devia- 
tions at stronger coupling. 
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B. Angle-resolved photoemission 

In photoemission experiments the intensity of photo- 
electrons is proportional to /(e)A(e,/c), where / is the 
Fermi distribution function and A is the spectral func- 
tion. In this section we discuss our results for the spec- 
tral function A{e^k) obtained from Eq. (33). For this 



via the relation 



purpose we follow a procedure commonly used to ana- 
lyze ARPES data: The superconducting state dispersion 
is compared to that of a reference dispersion in the nor- 
mal state (at a temperature well above Tc) and renor- 
malization effects are directly inferred by quantifying the 
differences. In addition, a linewidth analysis gives infor- 
mation about the imaginary parts of the diagonal self- 
energies. These procedures give naturally not as pre- 
cise information as an orbit ally resolved measurement. 
The reason is that the spectral function is obtained from 
the orbital-dependent self-energies via a matrix inversion 
[see Eqs. (29) and (33)], and all orbital contributions to 
the self-energy mix with each other when considering the 
spectral function. In other words, the self-energies are 
not diagonal in the band index, whereas the renormalized 
tight binding bands are not diagonal in the orbital index. 
This refers also to orbital-resolved photoemission experi- 
ments, as each orbital component of the spectral function 
is obtained by the same matrix inversion, Eq. (29), thus 



mixing orbital components of the self-energies. Neverthe- 
less, it is possible to compare theoretical and experimen- 
tal data when applying identical procedures to the theo- 
retically and experimentally obtained spectral functions. 
We call the related quantities effective self- energies. All 
results in this section are for g'^Xo — 1-1'^ x 10^/i|eV K. 



1. Effective self energies: dispersions and linewidth 

Motivated by experimental findings, we use a fitting 
procedure for our theoretically obtained spectral func- 
tions as discussed in the following. For fixed energy the 
momentum dependence of the spectral function (MDC: 
momentum distribution curve is peaked and often well 
approximated by a Lorentzian. In addition to MDCs, 
also the energy dependence for fixed momentum (EDC: 
energy distribution curve) can be useful to analyze the 
dispersion anomalies of superconductors that arise from 
coupling to bosonic modes. Whereas in the cuprates 
the EDCs exhibit a pronounced dip-hump structurej ^-^ * ^^ * 
there is no such pronounced behavior for iron-based 
superconductors.^ Here the self-energy effects are best 
extracted by examining the MDC-derived dispersions. 

When concentrating on such a Lorentzian peak cen- 
tered at momentum k and energy for a certain band 
well separated from all other bands, we consider the im- 
mediate region around this point in energy-momentum 
space. The energy dispersion of the considered band 
in relation to the bare band dispersion defines an effective 



(40) 



This equation is motivated by the analogous equation 
for a single band system, where describes the true 
low-energy part of the self-energy, and where the above 
equation represents an implicit equation for the renor- 
malized dispersion Sk- In the present case, we instead 
calculate the renormalized dispersion using the true self- 
energies, and then we define the effective self-energy 
via Eq. (|40|). 



By assuming that S'le is nearly constant in the mo- 
mentum region over which the Lorentzian describing the 
MDC spectral function is spread, the MDC spectral func- 
tion can be written in the form 



A{e,k) 



1 



n{i,{k-k,)r + {i:'l^{e,k)r' 
where ke is obtained from e = = Ske/^UE and 



(41) 



{^k - ^)\k^k. 



^LE — ^Le/^HE, 

dsk 
Ilk 



(42) 

{k-k,) = v,{k-ke), (43) 



and the group velocity is related to the effective self- 
energy by 



1 \ dCk/dk + dE'i^E/dk \ 



Zhe I i-d^ijdsk 



(44) 



The bare dispersion (k cannot be measured. However, 
it is temperature independent, such that the quantity 



S'(e) 



^k. 



(45) 



gives a measure of the self-energy effects only, S'(e) = 
^LE(^'^e)^ — I]le(^' ^e)^"""^- It cau be extracted experi- 
mentally by taking the difference between two MDC dis- 
persions at fixed k^. Here T^ef is a reference temperature 
that should be chosen deep in the normal state. 

Another interesting quantity, which determines the 
quasiparticle lifetime and is accessible by Eq. (41), is 
the linewidth function S''(e) = ^{^^(e, /Ce)/ZHE5 which 
we obtain by determining the full width at half maximum 
(FWHM) of the Lorentzian in reciprocal space, \Ske\^ and 
the group velocity in the direction of Sk^, 



t'\e) = v,3ke- 



(46) 



The procedure explained above is illustrated in Fig.[To] 
We choose T^ef = 50 K. Because the exact value of Zhe 
is unknown, we fix it at this temperature, = 1, 

and determine Zhe in the superconducting state so that 
the dispersions merge with the normal state reference at 
high energies. In Fig.[lO]the reference dispersion is shown 
along a cut in the first Brillouin zone for = 0. In the 
extracts (a) and (b) we present the dispersion branches 
of the ai^2 and the f3i bands. Below the Fermi energy 
we fitted the MDC dispersion to the spectral functions, 
which are shown by the black solid dots, and calculated 
the linewidth by the procedure explained above. 
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FIG. 10: Effective Self-energies: Illustration of the procedure 
used to extract self-energy effects from the electronic spectral 
function A(e, k) at Tref = 50 K. The width of A{e, k) for fixed 
energy determines the scattering rate ^ whereas the dif- 
ference of two dispersion relations at different temperatures 
gives the real part of the effective self-energy, S^ 
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FIG. 11: (a) Temperature dependence of the imaginary part 
Tj" of the effective self-energy at the ai band extracted from 
the cut in Fig. 10 'a). The black solid curve corresponds to 
the normal state at T — 50 K. The dashed blue, green, and 
red curves correspond to T = 13, 24, 31K respectively. The 
difference of the imaginary parts of the effective self-energy, 
S''(T)-S''(50K), at T = 13, 24, 31 K is shown at the bottom, 
(b) The same as in (a) along the cut at the /3i band which is 
depicted in Fig. [To|b). 



2. Quasiparticle lifetime in the normal state 

The imaginary part of the effective self-energy exhibits 
a linear dependence at high energies, 5]''(e) oc e_(Fig. 11), 
consistent with marginal Fermi liquid theory.!^ This re- 
sults from the coupling to the spin-fluctuation contin- 
uum, in particular from the slow decay towards high 
energies. The linear in energy contribution increases in 
magnitude with decreasing temperature, consistent with 
the temperature dependence of the spin-fluctuation con- 
tinuum. It persists at high binding energies in the super- 
conducting state. As we will see in the following section, 
in the superconducting state the appearance of the spin 
resonance leads to an additional feature at low energy, a 
small bump in the imaginary part of the self-energy seen 
in Fig. [TT] for T = 24K and, more pronounced, for T = 
13K. 



3. The superconducting state and the "kink^^ feature 

Upon entering the superconducting state the low- 
energy resonance in the dynamic magnetic susceptibil- 
ity appears below the particle-hole continuum. In this 
section we want to extract the effect of this resonance 
feature on the electronic structure. For this we compare 
the normal state and superconducting state dispersions 
and extract both, the real part as well as the imaginary 



part of the effective self-energy from those. The influ- 
ence of the resonance on the dispersion relation can be 
quantified using, e.g., the relation Eq. (45), 



where we obtain the renormalized dispersion relation 



;rT<Tc 



by fitting the Lorentzian in Eq. (41) to the theo- 



retically obtained spectral function in Eq. (33). That is 



a good approximation even in the superconducting state 
for energies not too close to the superconducting gap. 

Now the electronic dispersions = s^/Z^^ at the 
ai 2 bands and the /3i band can be determined. In 



Fig. 12 the temperature dependence of the real part of 
the effective self-energy is shown as obtained in Eq. (45). 



The appearance of the bosonic resonance leads to an 
effect on the electronic dispersion which is character- 
ized by the development of a peak in the real part of 
the effective self-energy, S'(e) (see Fig. 12), as well as 
a hump feature in the imaginary part of the effective 
self-energy T,"{e) = T,"{e^ke) (Fig. 11), which we see 
even more pronounced by having a look at the difference 
T."{T) - S''(Tref) [Fig.llija) and[lTJb) at the bottom]. 
Both effects follow the temperature dependence of the 
bosonic resonance. The coupling feature is situated at 
energies Aq + ^?^s5 where Aq is the gap at the corre- 
sponding Fermi surfaces that are approximately (i.e. on 
the scale ^~^) nested by an antiferromagnetic wave vec- 
tor Q. The broad maximum in Tj' which moves around 
^ 50 meV at all temperature is the result of coupling to 
the particle-hole continuum. 
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FIG. 12: (a) and (b) The temperature dependence of the 
real part of the effective self-energy extracted from the 
dispersion cuts shown in Figs. 10 'a) and|lQ[b). 



The fact that the above-discussed kink feature results 
from coupling to the spin-fluctuation resonance was al- 
ready predicted in ARPES experiments by Richard et 
al. There the authors show that the coupling feature 
follows an order-parameter-like temperature dependence 
which is consistent with the temperature dependence 
of the neutron resonance observed by inelastic neutron 
scattering. By comparing the electronic dispersion in 
the superconducting state with the one in the normal 
state (see the illustration in Fig. 10) they clearly show 
that the low-energy kink is emerging in the supercon- 
ducting phase. Within our calculations we are able to 
reproduce the shape and the absolute value of the real 
and imaginary part of the effective self-energy with a cou- 
pling strength of ^^xo = 1-17 x 10^/i|eVK. 

On the other hand, Koitzsch et al^ performed a de- 
tailed analysis of the temperature dependence of the elec- 
tronic dispersion relation deep into the normal state. 
The authors stress that there exists a bump feature 
at ~25 meV in the normal state which does not mod- 
ify at temperatures below 100 K, continuously moves to 
higher energies at higher temperatures and finally dis- 
appears above 200 K. Although a kink between 25 and 
40 meV would be consistent with the phonon spectrum^^ 
this would not explain the mentioned temperature de- 
pendence. However, there is a natural expanation of 
these observations within our model. Note that the 
spin-fluctuation continuum modeled by Eq. (pTj) and mo- 
tivated by neutron scattering experiments^^ exhibits a 
maximum at = ^^ax- This moves to higher ener- 
gies with increasing temperature whereas the amplitude 
of Eq. (11) continuously vanishes. According to that. 



the observations of Koitzsch et al. can be understood as 
the fingerprint of the spin fluctuation continuum, which 
persists in the normal state and weakens with increasing 
temperature. Energy scale, magnitude, and the temper- 
ature dependence are all consistent with this interpre- 
tation. The additional change of the kink feature when 
crossing Tc, as, e.g., clearly observed in the experiments 
by Richard et al.^^^ was not resolved in the experiments 
of Koitzsch et al. 



4. EDCs in the superconducting state 

Recent FLEX calculations have shown that the EDCs 
show a well-pronounced dip-hump feature.'^ This would 
essentially lead to a break in the dispersion relation, 
which is not the case for iron-based superconductors. In 
this section we want to take a look at the EDCs calcu- 
lated in our model. 

In Fig, p^a) and p^b) we compare the EDCs at dif- 
ferent cuts in the first Brillouin zone and at the two tem- 
peratures T = 15K and T = 31.0 K. As we have seen 
in the previous chapter, the less pronounced resonance 
feature at higher temperatures is also less pronounced 
in the effective self-energies (Fig. 11). Accordingly, we 



expect anomalous features, which can be related to the 
resonance, to appear at T = 15K but to be absent at 
T = 31 K. Actually, for the higher temperature the EDCs 
on the right-hand side of Fig. 13 'a) and IlsFb) show a 



smooth shape that should be compared to the left-hand 
side corresponding to the lower temperature. Here the 
electronic spectral function shows particular variations 
of this smooth lineshape, namely, small kinks (indicated 
by arrows). 

In Fig. [14] the superconducting spectral function along 
the F-X direction is shown. The maxima of the curves 



in Fig.[l3[a) and 13 ^b) give the EDC-derived dispersions 
indicated by the black triangles. We determine the su- 
perconducting gap at the pockets by the maxima of this 
dispersion, as illustrated in Fig. [T4j The superconduct- 
ing gap is nearly constant along the Fermi surface sheets 
with Ac,i ^ 11.6 meV, /^^^ ^ 8meV and /^^^ ^ 12meV. 

Another quantity of interest is the partial density of 
states (PDOS). Here the spectral function is integrated 
over a cut that crosses the Fermi surface, as we have 
illustrated in the inset of Fig. 15, i.e., PDOS(e) = 
Xl/c^G[- tt] ^'s(^) 1^2=0- There we present the PDOS 
curves at different ky values at the j3i pocket. For 




FIG. 13: Energy distribution curves at the ai,2 (a) as well 
the Pi pocket (b). We compare the EDCs at T = 15 K (left) 
and T = 31 K (right). The respective cuts are indicated in 
the inset. 
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FIG. 14: Superconducting state spectral function at T = 15 K 
along the T-X direction as well as the maxima of the EDO 
curves (black triangles) . The superconducting gaps are deter- 
mined as illustrated in the insets (A^ ~ 11.6 meV, ~ 
8meV, and A(3^ ^ ll.GmeV). 



which is absent at T = 31 K. The difference of this peak 
(at 25.5 meV) and the quasiparticle peak (here 12.6 meV) 
gives ~ 12.9 me V which corresponds to the resonance en- 
ergy nj^";^^^ ^ 12meV. 

In conclusion EDO's are not best suited to examine 
anomalies that arise due to the spin-fluctuation reso- 
nance because the features that can be referred to the 
mode are too small. However, they are essentially more 
pronounced in the partial density of states [PDOS(e)], 
as was already seen in experiment!^ and its derivative 

/ dPDOS(e) N 
V de )' 

C. Tunneling spectra and density of states 

Another useful tool to obtain information about the 
density of states (DOS) is scanning tunneling spec- 
troscopy (STS). 

We are interested in the spin conserving tunneling cur- 
rent when applying a voltage V. The tunneling current 
in this case is given by 



PDOS{e), ' r'= 15K 



I{eV) = \To\ 
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FIG. 15: Partial density of states, i.e., the added spectral 
function X^/.^^^ ^| Ak{e)\kz^o- We present the PDOS at dif- 
ferent /cy -values at the /Si band, as illustrated by the inset. 
The order of the PDOS curves from top to bottom corre- 
sponds to the order along the arrow in the inset. At T = 15 K 
the curves exhibit a dip at ^ 25 — 30 meV which is absent at 
T = 31K. 



T = 15 K they clearly exhibit a dip feature at an en- 
ergy of Qres + Aq ~ 25 — 30 meV, which is absent at the 
higher temperature T = 31 K. 

In order to enhance the coupling feature a common 
procedure is to take the derivative of the partial den- 



sity of states. This quantity is shown in Fig. [16] and we 
see that ^^^^^^^^ exhibits a clear maxima at T = 15 K 



de 
2^ 



M^)^2{e^eV) [/(e) -/(e + e^)], 

(47) 

where i^i,2(e) = Xlfc ^1,2(^7^) are the densities of 
states computed from the spectral functions ^41^2 of the 
two materials in contact, and / is the Fermi distribu- 
tion function. We assume here for simplicity an isotropic 
tunneling matrix element Tq. 

For a superconductor-insulator-normal metal (SIN) 
tunneling junction the energy dependence of the normal 
metal density of states in the vicinity of the Fermi en- 
ergy, Ci?, is negligible. In this case one defines the matrix 




dPDOS 








de 




f 
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FIG. 16: Partial density of states PDOS(e) curves as in Fig. 15 
at T = 15K (left) and at T = 31 K (right) and its derivative 
The dip in the PDOS leads to a peak in its deriva- 
tive which is absent at T = 31 K. 
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element |Mop = \Tq\^ vn{^f)i where vn is the local den- 
sity of states in the normal metal. The SIN tunneling 
conductance is then given by 



psiN^eV) 



= e\Mc 



dis 



IN 



dV 

oo 



(48) 



de_ us{^) 

27r 4/cBTcos2(el//2/cBT)' 



with the density of states ^^^(e) = ^^(e, /c), where 
As is the spectral function of the superconductor. For 
the coupling strength g'^xo = 1-17 x 10^ /i^eVi^, 
which is the value we obtain from our comparison with 
ARPES experiments, we present the tunneling conduc- 
tance Psin{^) in Fig- 17 for different temperatures. 

Van Hove singularities at specific points in the Bril- 
louin zone can cause an intensified scattering with spin 
excitations, leading, e.g., to a well-pronounced break fea- 
ture in the spectral function and therefore a suppression 
of the density of states i n the aff ected energy range in 
cuprate superconductors^^l^^I^^I^. This however is ab- 
sent in our calculations. 

Our spectra show a strong particle-hole asymmetry 
that has been also seen in experimental data.^^ 

Nevertheless, it is interesting to investigate how the 
tunneling conductance is modified by changes of the 
coupling strength and when significant features appear 
at low energies that can be attributed to the effect of 
spin fiuctuations. Therefore we calculated the differen- 
tial tunneling conductance dl/dV for different coupling 




FIG. 17: SIN tunneling conductance at different temperatures 
in the range 13 — 37i^ in the superconducting state for the 
coupling constant g^Xo — 1-17 x 1^^ ji^eV K. The curves 
show no clear hint of a coupling to spin excitations, such as a 
dip feature in the case of cuprate superconductors. 



strengths g'^xo = 0.68 ... 2.42 x 10^ /i| eVK for a SIN as 
well as a SIS junction. 

Both psiN and psis are presented in Fig. [18] The 
SIN tunneling spectra show the development of a dip- 
hump feature for increasing coupling strength, as in the 
case for cuprate superconductors. As can be seen from 
Fig. [Tsj for SIN tunneling the coupling predominantly 
influences the occupied side, whereas the unoccupied side 
remains nearly unaffected. In contrast, the SIS spectra 
are symmetric and therefore the dip feature appears on 
both sides. For a coupling strength of g'^Xo — 3.44 x 
10^/i|eVK we see a clear dip at ~ SOmeF that can 
be traced back to the resonance in the spin-fluctuation 




FIG. 18: SIN (top) and SIS (bottom) tunneling conduc- 
tance for different coupling constants g^x^ — 0.68. ..3. 44 x 



lO^li^eV K in units e|Mo|^ and e|To 
temperature T = 15 



respectively at the 
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spectrum. 



VII. ACKNOWLEDGEMENTS 



VI. SUMMARY AND CONCLUSIONS 

In this work we presented a model of spin fluctuations 
coupling to electrons in the normal as well as the su- 
perconducting state. We find that the high-energy spin 
fiuctuation continuum provides imporant contributions 
to the pairing interaction, leads to a linear in energy 
broadening of the quasiparticles, and contributes to the 
renormalization effects of the dispersion both in the nor- 
mal and in the superconducting state. In the normal 
state it leads to a broadened kink feature in the disper- 
sion, similar to that observed in experiments by Koitzsch 
et a/.^In the superconducting state, it leads to a sharper 
kink feature as experimentally observed by Richard et 
al^ and by Wray et al]^ 

The superconducting gap originating from the low- 
energy (<200 meV) part of the dynamic susceptibility 
supports an state, as proposed earlier.'^ The in- 
vestigation of the scattering rate of quasiparticles and 
spin fiuctuations shows that the renormalization effects 
are strongest at the approximately nested Fermi surface 
sheets and exhibit characteristics of these. 

The sharper kink feature at certain energies in the 
eletronic dispersion in the superconducting state is due 
to the low-energy resonance in the spin-excitation spec- 
trum. It can be quantified following experimental proce- 
dures to extract them directly from the spectral function. 
A comparison with the ARPES experiment by Richard 
et al^ enables us to estimate the strength of the cou- 
pling between electrons and spin fiuctuations. We find an 
intermediate value, certainly smaller than an analogous 
analysis in the cuprates yields. 

We discuss the renormalization of the Fermi surface 
and find that coupling to spin fiuctuations leads to a shift 
of the chemical potential to higher values by ~ 20 meV, 
accompanied by a shrinking of all Fermi surfaces. The 
relation between the chemical potential and the charge 
carrier density is linear for small to moderate values, and 
non-linear for higher coupling strengths. 

Investigating the energy dependence of the spectral 
function as well as the partial density of states and the 
differential tunneling conductance we find that for the 
coupling strength necessary to reproduce the self-energies 
of the ARPES experiments, there are no such pronounced 
coupling features in the tunneling spectra as observed in 
the case of cuprate superconductors. However, in the 
partial density of states, obtained when integrating the 
spectral functions perpendicular to the Fermi surface, 
coupling features are clearly visible. 

We conclude that the coupling between spin fiuctu- 
ations and electronic excitations provides an important 
mechanism for superconductivity in iron pnictides. 



RG would like to thank the Karlsruhe House of Young 
Scientists for financial support during his stay at Royal 
Holloway, University of London. The work of AH was 
generously supported by the South East Physics network 
(SEP net) during his stay at Royal Holloway, University 
of London. 



Appendix A: Energy and momentum dependence of 
the diagonal and ofF-diagonal self-energies 

In this appendix we discuss the diagonal self-energies 
entering Eq. (30) and Eq. (32). The energy and mo- 



mentum dependence of the renormalization function Z is 
presented in Fig. 19 At low energies, the renormaliza- 



tion function takes values of about 2 {yz channel) and 
4 {xy channel), and at higher energies it takes values of 
about 1.5 {yz channel) and 2 {xy channel). Results for 
the xz channel are analogous to the results for the yz 
channel. 



Also shown in Fig. 19 is the particle- hole asymmet- 
ric part of the diagonal self-energy, corresponding to 
Eq. (30). Looking at the xy channel of the real part 
we see a drastic change for energies |e| > \ fires + ^7! at 
{kx^ ky, kz) ^ (0, ±7r,0). This corresponds to the ap- 
pearance of the shallow electron pocket at the {3 bands. 



In Figs. |20] and [21] the energy dependence of the real 
and imaginary parts of the off-diagonal self-energy as 
well as the renormalized gap is shown. The gap func- 



tion, shown in Fig. 21, is real in the energy range <15 
meV, given by the energy of the magnetic resonance. The 
gap itself is of comparable magnitude. Above this energy 
range, the gap function acquires a considerable imaginary 
part, which is of the same order of magnitude as the order 
parameter itself. Note that these effects are very small in 
the yz channel near k = (±7r,0, /c^), i.e., at the electron 
pockets, and in the xy channel near k = (0,0, /c^), i.e. 
at the large hole pockets. Thus, strong coupling features 
in the gap function are almost purely in the xy channel 
for the electron pockets and in the xz and yz channel for 
the large hole pockets. The resonance imprints at ener- 
gies fires + and fires + ^7 also are indicated in Fig. [6] 
i^Ies^^^ ~ 12meV). This effect is most pronounced in 
however, due to the enhancement of the renormal- 



ization function Zn at these energies (see Fig. 19) the 



renormalized gap weakly varies, as is especially seen 
in the yz channel in the upper panel of Fig. [21] 



Appendix B: Local Sum Rule 

In this Appendix we shortly review a sum rule 
that holds for the two point correlation function 
{S^{t,Ri)S^{t',Rj)), where S^{t,Ri) denotes the spa- 
tial a = x^ y^ z component of the spin at site Ri and in 
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FIG. 19: Energy and momentum dependence of the real part 
of the renormahzation function Zn and of t he real part of the 
dispersion correction according to Eq. (30), (C^ — C,^)/Zn = 
(T^+T^) /2Zn, at T = 15 K. The energy dependence is shown 
along the momentum-cut {ky G [— tt : tt), /c^ = /c^ = 0}. 



orbital m. Then the spin structure factor is defined by 



J — ( 

EE 



dt e 
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FIG. 20: Energy and momentum dependence of the real and 
imaginary part of the off-diagonal self-energy $n at T = 15 K 
along the cut {kx G [—tt : tt), ky = kz = 0}. 



spins. This is realized in the Heisenberg limit, where 
the hopping of electrons between the different sites is 
suppressed. However, in more itinerant systems (metallic 
regime) correlations between the spins are reduced and 
the integrated weight should be smaller. However it was 
proven in experiment that the sum rule also holds in the 
case of iron pnictides.^^ 



Integrating this quantity over frequency yields the 
static spin structure factor 

^ mn 

Further integrating out the momentum and taking the 
trace gives 

E / 7S^^""(^) = ^EE(^™(^)^"(^^)) 
= ^E(^(^')-^(^')) 

i 

= S{S^1). 

Here we defined the spin S{Ri) = Sm{Ri) on each 
site i. The last equality only holds for a system of equal 



Appendix C: High energy renormalization factor 



Starting from Eq. (41) we have introduced a high- 



energy renormalization factor Zhe = 1 — ^he(^)/^ 
order to treat the renormalization originating from the 
high-energy incoherent part of the spin fluctuation spec- 
trum. In the procedure explained above we obtained a 
linear temperature dependence Zhe = a -\- ^T. In this 
section we provide a motivation for such a temperature 
dependence. 

By introducing a high-energy cutoff Uc in the bosonic 
spectrum, the self-energies consist of two contributions, 
E^(e,it) = ^lf{e,k) + E™(e,it). Whereas the low- 
energy contribution can be calculated exactly, we have 
to treat phenomenologically. As we have seen in 

Sec. |IIB[ the sum rule for the susceptibility implies 
that f_^duj J d^qS{uj^ q) remains constant and tempera- 
ture independent. As we assumed in our model identical 
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FIG. 21: Energy and momentum dependence of the real and 
imaginary part of the renormaUzed order parameter = 
(A^ + $^)/Zn at T = 15 K along the cut {A:^ G [-tt : tt), ky = 
= 0}. 



weights for the orbital contributions Xn{^^k)^ the sum 
rule applies to each of them separately. Introducing the 
high-energy cutoff we obtain 



/ 



yJ\^\<uj, l-exp(-w/T) 

Imx„(w,q) 



/ 



i\oj\>u. 
^LE + ^HE = I 



1 -exp(-a;/T) / 
const. (CI) 



Whereas the low energy behavior of the susceptibility 
is motivated by recent experimentJ^, the high energy 
spectrum is unknown. On phenomenological grounds we 
expect it to be broadened in momentum and to decay 
at high energies. Furthermore the temperature under 
consideration is small compared to the energy cutoff, i.e. 
uOc/T » 1 and therefore 



I- I (fq 



d^ , ^"^^"("'^L (C2) 



\uj\<UJc 



1 - exp(-cj/T) 



By numerically integrating the second term on the right- 
hand side we find a linear dependence in temperature, 
^LE = ^LE - 7le^ (see Fig. [22| implying that 
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FIG. 22: Value of the low-energy integrated weight /le = 
^LE — TleT" as function of temperature T, showing a linear 
dependence. 



f 



dwlmx„(w) / - 7?LE + 1-LE.T- 



(C3) 



Let us now consider the high-energy contribution of the 
self-energy. We start from an expression for the electronic 
self-energ}i^ of electrons coupling to the spin fluctuation 
mode (in real gauge). 



vR/ , N _ o^2Y^^ ImG^^(C, k')lmxn{^^ k - k') 

k' ojX 



tanh ( ^ ) +coth 

\2T ) \2TJ 



do; 
oo 27r ' 



with the abbreviations Ylk' = I = L 

and = We are interested in the high-energy 

contribution of the spin susceptibility, for \uj\ > c^c and 
therefore restrict the integration over u in the expression 
above to this range. At high energies the spin- fluctuation 
contribution will become more and more local, so that it 
is well approximated by a momentum eigenfunction ex- 
pansion with a few eigenfunctions r]s{k) resulting from 
close neighbors. At the same time, we relax the ap- 
proximation that the spin susceptibility only depends on 
q = k — k' . Thus, we write 

Xn{^,k- k') Xn{^,k,k') = ^r]s{k)xn,s{^)Vs{k'y 

s 

with Xn,s(^) = Ylk,k'Vs{kyXn{^,k,k')r]s{k')' The 
eigenfunctions can be classified according to the ir- 
reducible representations of the point group, and 
they are assumed to be orthogonal and normalised, 
Ylk'^s{kyr]s'{k) = Sss'^ and built to a complete set, 
^gT]s{k)r]s{k^y = ^k,k'' Thus, the corresponding com- 
ponents of the self-energy, E^g(e) = X^jr. ^s(^)*S^(e, /c), 
are given by 
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^HE 



tanh 



2T 



coth(|,) 



with Gl,,(C) = E'k'Vsik'rG^'iCk'). With the 
partial density of states i^n,s(C) = ~:^I^^nn s(C) we ob- 
tain 



The second line in this expression ensures that |cj + CI is 
always of order cjc, as it mostly contributes for ( > —2T 
for positive uj and ( < 2T for negative uj. As ^ T, 
for |e| <^ uJc this allows us to expand the denominator in 
the first line in e, leading to 



, dcj r 



^n,8(C)ImXn,8(^) 

e - (cj + C) + ^0 



tanh 



2T 



- sign(a;) 



c). (C4) 



tanh 



with i/n,s±(C) = [i^n,s(C)±^n,s(-C)]/2, and whcrc wc uscd 
ImXn,s(— ^) = — IniXn,s(^)- The ^-integral gives up to 
second order in a power expansion in temperature 



-(C) 



(C)e 



oj + C 



,(0) 



{uj - 



.(0)6 



with a numerical constant a of order 1. We see that the 
temperature-dependent terms are small compared to the 
temperature-independent terms. The terms independent 
of e contribute to the band renormalization. If all orbitals 
would contribute equal, it would correspond to a constant 
shift of the bands. The leading contribution is 



/■HE 



-45^ 



duj 
2^ 



Imxn,s(^) / dC 



-(C) 



.(C5) 



For different orbital contributions to the spin suscep- 
tibility, this quantity would lead to orbital-dependent 
band shifts, and thus to a slightly temperature-dependent 
renormalization of bands. The temperature-independent 
part can be absorbed in the definition of the tight binding 
bands and the orbital as well as Brillouin zone average 
in the temperature-dependent chemical potential. The 
remaining contributions are neglected in this paper. 

The terms linear in energy contribute to the high- 
energy renormalization factor Z^^. We can write in lead- 
ing order 



7HE 



^ J ^^"^^^'^(^) ^C 



^n,.+ (C)^^^g^ 



For our calculation we restrict with regard to the di- 
agonal self-energy contributions to the Brillouin zone av- 
erages (corresponding to the basis function 5 = with 
r]o{k) = 1), and we define the orbital average Zhe = 
(Z^q)^. According to that the temperature dependence 
of the high-energy contribution Zhe is determined by the 
temperature dependence of XhE' leading to 



Appendix D: Kramers-Kronig Relations 



(C7) 



For the diagonal self-energies in Eq. (24) 



we obtain the following Kramers-Kronig relation, 
Re[E^(e,fc)-E^(oo,fc)] = -V T de'^^^^^^. 



which implies. 



1 

Rep^(0,/c)-S^(oo,/c)] = -V / de' 

J-00 



, ImS^(e^/c) 



Here, V denotes the principle value integral, and 
11^(00, k) = lim|g|^oo (^5 ^)- With the definition 

S±(e,fc) = [S^(e,fc)±S^(-e,fc)]/2 
we can write down corresponding relations for the renor- 
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malization function and band renormalization, 

, ImS+(e',fe) 



ReZ„{e,k) = 1 



Jo 



de' 



e'2 - e2 



2e2„ de' ImS-{e',fe) 



+ 



Jo e 



Note that Iml]+(e, k) < for all e and k (whereas ImS^ 
can have either sign), and consequently 

ReZr^iO^k) = 1-- Hde^ ImS+y,fc) 
Jo 

2 /""^ 

lim ReZ^(e,fc) = 1 + ^ / de' Imi;+(e', fc) < 1, 

\e\^oo Tre^ Jo 

Finally, for completeness we also present relations for the 



band renormalization, 

ReC;(0,fe) = CM(fe) + ReS^(0,fe) 

lim ReC;(e,fc) = ^k) + Re^^iO,k) 

|e|-!-oo 



2 r 



de 



, Iml]-(e',fe) 



In our case, ImE~ is predominantly positive, leading to 
negative corrections at large energies. Finally, we have 
the sum rule 



f 

J — c 



de'e'A{e',fc) = ^ |a;:{fe)P lim ReC;:(e,fe) 
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